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We numerically study controllability of quantum dynamics in a perturbed one-dimensional double- 
well potential by using an optimal control theory. When the perturbation strength is small the 
dynamics of an initially localized Gaussian wavepacket shows coherent oscillation between the wells. 
It is found that as the increase of the strength and/or the number of frequency components of 
the perturbation the coherent motion of the Gaussian wavepacket changes to an irregular one with 
irreversible delocalization. We investigate the controllability of the system depending on the per- 
turbation parameters and the initial quantum state by focusing mainly on the delocalized state 
generated by the polychromatical perturbation. In relatively long-time control for the Gaussian 
wavepacket and the delocalized state, we show that it is well-controllable via the first excited state 
doublet in spite of the perturbation parameters. On the other hand, in the relatively short-time 
control we show difficulty of the control for the delocalized state because of the numerous local 
minima. Furthermore, it is demonstrated that chaotic phase space structure of the system assists 
the controllability of the delocalized state. 

PACS numbers: 05.45.M+, 05.45.Gg, 66.35.-l-a 



I. INTRODUCTION 



Decoherence of the quantum dynamics and a control of 
the wavepacket dynamics by using an external force field 
are recent attractive topics. There are many studies on 
environment-induced quantum decoherence by coupling 
the quantum system to a reservoir 

0,11 nil- Quantum 
dissipation due to the interaction with chaotic degrees 
of freedom has been also studied |E IE • It is expected 
that the quantum dynamics is sensitive to any disruption 
of coherence as it occurs due to an unavoidable coupling 
to the environment. One of obstacles for experimental 
realization of a quantum computer is decoherence of a 
quantum state caused by coupling to the other degrees 
of freedom 0. 

On the other hand, in the experimental realization of 
control of the chemical reactions and the optical molec- 
ular devices, a design of the control field based on op- 
timal control theory (OCT), is also an important topic 
miilllllllll- Recently, some methods to control the 
transition between quantum states have been proposed 
CllllllllllllllilllSIIllll. Our main interest is in 
the controllability between the quantum states in, what 
we call, a quantum chaos system 24] . 



the following Hamiltonian, 



In the previous papers |25|, |2g , we investigated influ- 
ence of chaos on quantum tunneling and decoherence in 
a parametrically driven double-well system modeled by 
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where q and p represent the generalized coordinate and 
the conjugate momentum, respectively. Here M denotes 
the number of mutually incommensurate frequency com- 
ponents {riiji = 1,2, M}. We choose off-resonant fre- 
quencies which are far from both classical and quantum 
resonance in the corresponding unperturbed case (e = 0). 
The number M of frequency components in the poly- 
chromatic oscillation corresponds to the number of de- 
grees of freedom coupled to the unperturbed double-well 
system ^7J. The details of relation between the time- 
dependent Hamiltonian and the autonomous representa- 
tion have been shown in Ref . ^3 ■ 

Note that in the previous works of Lin et ai, they 
dealt with a double-well system driven by single forced 
oscillator, therefore, the asymmetry of the potential play 
a role in the chaotic behavior and tunneling transition 
between the quasienergy states |2E|2E|33- However, in 
our model the potential form is remained symmetric dur- 
ing the time evolution process, and different mechanism 
from the forced oscillation makes the classical chaotic be- 
havior. We classified the motions of an initially localized 
Gaussian wavepacket depending on the perturbation pa- 
rameters (M, e), i.e. coherent motions like an instanton 
and irregular motions of a delocalized state 26]. We 
call the delocalized state in the irregular motion chaos- 
induced delocalized state in this paper. 
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The purpose of this article is to present numerical re- 
sults on the control of the quantum dynamics in the one- 
dimensional double-well potential with the coherent per- 
turbation. The controllability of the system depends on 
the perturbation parameters and the initial state. In the 
relatively long-time control the displacement of the Gaus- 
sian wavepacket is well-controllable by a designed field 
based on OCT in the perturbed double-well potential re- 
gardless of the external parameters. Furthermore, in the 
relatively short-time control we show the difficulty of the 
control for the chaos-induced delocalized state generated 
under the perturbation because of numerous local opti- 
mal fields, compared to the case of Gaussian wavepacket. 
It is numerically demonstrated that when the system is 
fully chaotic by the increase of the perturbation strength 
e the controllability of the delocalized state can be en- 
hanced by the chaotic behavior. We call the phenomenon 
chaos-assisted control in this paper. 

The remained sections of the present paper are as fol- 
lows. In the next section we give some backgrounds and 
the numerical setting of the model and brief review of 
chaos-induced delocalized state. In Sect. 3 we shortly 
explain the optimal control theory used in the numeri- 
cal calculation. In Sect. 4 we give numerical results for 
control between Gaussian wavepackets in the polychro- 
matically perturbed system. In Sect. 5 we show results of 
the transition between a chaos-induced delocalized state 
and a Gaussian wavepacket. In Sect. 6 and 7 the dif- 
ficulty of the control due to the local minima and the 
chaos-assisted control are investigated for the short-time 
control, respectively. The last section contains summary 
and discussion. 



II. MODEL AND CHAOS-INDUCED 
DELOCALIZED STATE 

In this section we give the details of the setting in the 
Eq.(2) and in the numerical calculation. Furthermore, 
we give a brief review on occurrence of a irregular mo- 
tion in the delocalized state. The parameters are set as 
a = 5, e = 0.1 ~ 1.0, to emphasize the tunneling effect 
in the energetically and/or dynamically forbidden region 
during the time evolution. We used the following Gaus- 
sian wavepacket localized at a bottom of the right well 
as, 

^.(g) = (a7r)i/4exp{-^i^^}, (3) 

where go = is the postion of the right bottom and 
the spread of the initial packet a ~ 0.3 as the initial 
state at t — 0. h — 1.0 and time mesh St is order 
of 10~^. The spatial discretization is chosen so that 
2®(= 256) points cover the interval [—5.5,5.5]. The 
Gaussian wavepacket can be approximately generated by 
the linear combination of the ground state doublet as 
V'j('7) ~ "72 ('/'o (9) + 4'i{'i))i where and (pi denote the 



ground state doublet. (See Fig.l.) Indeed, the ammo- 
nia molecule is well described by this barrier height with 
two doublets with energies below the top of the barrier, 
and the quantum dynamics is well approximated by the 
several levels from the bottom of the well in the weekly 
perturbed case [sH- 
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FIG. 1: The initial (target) Gaussian wave function \4>i{t = 
0)P {\'<Pf{t = tf)? = \^f? ), respectively. The dotted line 
shows a delocalized state as the initial state, which is gener- 
ated by time-evolution {t = 5.54 x 10^) in the perturbed case 
[M — 10, e = 1.0). The curve of the double-well potential is 
also plotted. 

In the classical dynamics, such a system shows chaotic 
behavior by the oscillatory force A(t) ,25|,i22|. The New- 
ton's equation of the motion in the monochromatically 
perturbed case {M = 1) is 

(fa 

^ - (a - e sin nt)q + = 0. (4) 

Note that the classical system is known as nonlinear 
Mathieu equation which can be derived from surface 
acoustic wave in piezoelectric solid j32j and nanomechan- 
ical amplifier in micronscale devices |33l |. In the poly- 
chromatically perturbed cases (M > 1) the smaller the 
strength e can generate chaotic behavior of the classical 
trajectories the larger M is |!26l|. 

We define transition probability PL{t) of finding the 
wave packet in the left well and a degree of coherent 
motion, AP^, based on the fluctuation of the transition 
probability as, 

PL{t) = f |V'(g,i)prfg, (5) 

J —00 

APi ^ ^{{PL{t)-{PL{t))Tf)T. (6) 

where (...)t represents time average value for a period 
T(~ 9.4 X 10^) liil. In the cases that perturbation 
strength e is relatively small, Phit) can be interpreted 
as the tunneling probability that the initially localized 
wave packet goes through the central energy barrier and 
reaches the left well. 
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FIG. 2: (a)Transition probability PL(t) as a function of time 
t for M — 10 with e = 0.1, 1.0. (b)e dependence of the degree 
of the coherent motion APl in the case. 



Figure 2(a) shows the time-dependence of PL{t) in the 
polychromatically perturbed case (M = 10). Apparently 
we can observe two kinds of time-dependence, i.e. coher- 
ent (e = 0.1) and irregular (e — 1.0) motions. In order 
to estimate quantitatively the difference between coher- 
ent and irregular motions, we use the degree of coherent 
motion AP/, in Eq.(6) and roughly divide the types of 
motion of wavepacket as seen in Fig.2(b). In the coherent 
motions, the values of AP^'s are almost same to the un- 
perturbed case, i.e. APj^ > 0.3. In the irregular motions 
which are similar to the stochastically perturbed case, the 
values of AP^'s becomes much smaller, i.e. AP^ < 0.2. 
In the irregular motions, the wavepacket is delocalized 
over the both the wells as seen in Fig.l. The delocal- 
ized state cannot show the recurrence to the Gaussian 
shape within the accessible computation time. We call 
delocalized states in the irregular motions chaos-induced 
delocalized state in this paper. 

It is well-known that an existence of a turning point 
breaks the Gaussian shape during the tunneling process 
even in a coherent motion. However, in the coherent mo- 
tion after the tunneling process the wavepacket goes back 
to the Gaussian shape again. In the following sections we 
give numerical results of the controllability for the Gaus- 
sian wavepacket and the chaos-induced delocalized state 
corresponding to the irregular motion. 



system, i.e. e = 0, as the controlled system the numeri- 
cal result does not change in the essential point, except 
for one in Sect. VII where the effect of the perturbation 
strength is discussed. 

The OCT gives an optimized external field E{t) in or- 
der to produce the target state $/ at the target time 
t = tf starting from an initial state ipilt = 0). The ob- 
jective functional J[E{t)] so as to maximize is 



J[E{t)]^\{Mtf)\^f)\-c^ / Eitydt 

Jo 



2Re 



{Mtf)\^f) l'\i^f{t)\-l-^-iHit)\ut))dt 



Here 'ipi{t) is a wave function driven by the optimal field 
E{t) which has been determined by the quantum states, 
ijjiit) and ipf{t) as. 



(9) 



Tpfit) is an inversely evolving quantum state starting 
from the target state ■(/'/(t = tf) = (f>/, which can be in- 
terpreted as a Lagrange multiplier to satisfy Schrodinger 
equation. The second term in Eq.(8) is put to limit the 
resulting field intensity, where a is a penalty factor to 
weight the significance of the external field. The optimal 
field strongly depends on the initial and final states, the 
potential parameters and the target time tf. Note that 
the existence of the optimal field has been proven for 
large target time and small a limit in the unperturbed 
case 1^. Generally, when we can optimize the transition 
between the energy eigenstates by a tt— pulse, the tar- 
get time must be at least i/ ~ 0{h/AE) to resolve the 
energy difference AE. 

The relation between the controllability and chaotic 
behavior in quantum systems has been investigated in a 
context of OCT H. 



III. OPTIMAL CONTROL OF QUANTUM 
DYNAMICS 

In this section, we briefly explain the OCT |^ we 
used in the numerical calculation. 

We consider the total Hamiltonian H{t) combined the 
system -ffo(i) with the interaction by the external field 
E{t) as, 

H{t) = H,{t)-^i{q)E{t). (7) 

The second term in Eq.(7) means the interaction between 
the transition dipole moment /i(g)(= q) and the exter- 
nal field E{t), where the charge is set unity. Note that 
we adapted the Hamiltonian Ho(t), including the time- 

2 

dependent part —A{t)^, as a controlled system. How- 
ever, we confirmed that if we use the static double-well 



IV. CONTROLLABILITY OF THE GAUSSIAN 
WAVEPACKET 

In this section, we give numerical results that displace- 
ment of the Gaussian wavepacket can be well controlled 
in the perturbed double-well potential. 



A. long-time control {tf — 300) 

We would like to control the transition between non- 
stationary states. Here we take Gaussian wavepackets 
localized around q = qaiq = —qo) as the initial(target) 
state, respectively. (See Fig.l.) Accordingly, in the 
quantum control there is no energy transfer between 
the initial and target states. We take the target time 
tf — 300 that is much smaller than a tunneling time 
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FIG. 3: Long-time optimal control of the Gaussian 
wavepacket in the monochromatically perturbed case{M = 
l,e = 0.1). (a)Time-dependence of the squared final over- 
lap Jo{t) after 20 iterations, (b) Convergence of the optimized 
objective functional values J and Jo as a function of the num- 
ber of iteration. The panels (c) and (d) show the optimized 
fields E{t) and the power spectrum I{uj) of the field. We used 
a = 1.0 and tf — 300 in the calculation. 



scale T = h/AEi2 ^ 9.7 x 10^ in the unperturbed case 
(e = 0). 

Figure 3(a) shows the final overlap Jo(i) = 
|(f/'i(i)|$/)P after 20 iterations in the monochromatically 
perturbed case (M — 1) for the target time tf = 300. 
As we can expect, based on the OCT, the iteration of 
the Jo shows fast convergence [? ]. The final state by 
the optimal field is well overlapped for the target Gaus- 
sian function. In Fig. 3(c) and (d), time dependence of 
the optimized field and the power spectrum are shown. 
The main peak of the power spectrum corresponds to 
energy difference between the ground state and first ex- 
cited state doublets that energy is lower than the barrier 
height. We can immediately confirm that the wavepacket 
travels between the wells via the first excited state dou- 
blet by calculating some quantum mechanical mean val- 
ues, < H >,< q >,<p >, and so on during the time- 
evolution process. 

Figure 4 shows the relatively long-time optimal control 
for the displacement of the Gaussian wavepacket in the 
polychromatically perturbed cases {M = 10). The ftmc- 
tion form of the optimal field depends on the extent of the 
chaotic behavior in the classical system. In the strongly 
chaotic case (e — 1.0) where the perturbation strength 
is large, the optimized field becomes a complicated func- 
tion of time with many modes. The peak structure of the 
power spectrum shows that a lot of modes are necessary 
to shift the Gaussian wave packet in the strongly chaotic 
cases. 

As a result, for the relatively long-time control between 
the Gaussian wavepackets with no energy deference it is 
almost perfectly controlled (Jo(t/) ~ 1.0) regardless of 
the perturbation type. 



FIG. 4: Long-time optimal control of the Gaussian 
wavepacket in the polychromatically perturbed case (M = 
10). (a) and (c) are the optimal fields after 20 iterations for 
e = 0.1 and e = 1.0, respectively, (b) and (d) are the power 
spectra for e — 0.1 and e — 1.0, respectively. We used a — 1.0 
and tf — 300 in the calculation. Jo(t/) ~ 1.0 in the both 
cases. 



B. short-time control (tf = 30) 

In this subsection we also consider the control prob- 
lem between the Gaussian wavepackets in the relatively 
short-target-time tf = 30. It can be expected that the 
optimal field has much different function form from the 
TT— pulse type for the short target time tf, and that the 
local minima may cause instability of the convergence to 
the global optimal field. 
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FIG. 5: Short-time optimal control of the Gaussian 
wavepacket in the monochromatically perturbed case {M — 
l,e = 0.1). (a)Time-dependence of the squared final overlap 
Jo{t) after 100 iterations, (b) Convergence of the optimized 
objective functional values J and Jo as a function of the num- 
ber of iteration. The panels (c) and (d) show the optimized 
field E{t) and the power spectrum I{uj) of the field, respec- 
tively. We used a = 0.1 and tf = 30 in the calculations. 
Jo{tf) ~ 0.9 after 100 iterations. 
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FIG. 6: (a) The probability density of the target Gaussian 
wavepacket and the final wavepacket achieved by using the 
optimal field in Fig. 5(c). (b) The time-dependence of the 
uncertainty product ApAg during the time-evolution by the 
optimal field. 



In Fig. 5 we show the result of the short-time optimal 
control in the monochromatically perturbed case {M = 
1, e = 0.1). By adjusting the penalty parameter a for the 
optimization the optimal field can be almost found out 
even for the short-time control (Jo(t/) ~ 0.9). 

However, as seen in Fig. 6(a) the final wavepacket has 
uncontrollable tail even after many iterations. The short- 
time optimal control is much harder than the relatively 
long-time ones. Figure 6(b) shows the time-dependence 
of the uncertainty product ApAq = y/{{q— {(i)Y)\/{{p — 
(p))^), during the time evolution by the optimal field in 
Fig. 5(c), where (...) denotes quantum mechanical aver- 
age. We use uncertainty product as a simple measure 
for quantum fluctuation in the phase space. It is found 
that the Gaussian wavepacket does not delocalize dur- 
ing the process. Furthermore, we have confirmed that 
the similar optimal short-time control can be achieved 
for the other cases with different perturbation parame- 
ters (M = 10,e = 0.1), (M = 10,e = 1.0), and so on. 
Here, we used a = 0.1 to adjust the strength of the ex- 
ternal field E{t). The small value of the penalty factor 
allows the optimized filed with the large amplitude. The 
existence of a lot of different optimized fields depending 
on a means the existence of a lot of local minima in the 
solution space. The details of the a— dependence will be 
given elsewhere [46l |. 

Generally speaking, the non-stationary state ^'(t) at 
time t can be expressed by a set of the expansion co- 
efficients {cj}t as '^{t) — '^jCj{t)(l)j, where are 
eigen states in the unperturbed case. Accordingly, the 
quantum control from the initial state at t = to the 
final state aX t = tf is a, problem to answer a question 
that how we can get the set of the coefficients {cj}t=t} 
in $^ = '^jCj{t = tf)(j)j in starting from {cj}t=o with 
obeying the Schrodinger equation. We can see that the 
transition between the Gaussian wavepackets can be al- 
most attained by a change of the relatively simple phase- 
relation between the states via the first excited state 
doublet as (co(0) = 1/^/2, ci(0) = 1/^/2) ^ (co(t/) = 
l/y2,ci(i^0 = -l/V2). 



V. CONTROLLABILITY OF THE 
DELOCALIZED STATE 

In this section, we investigate the controllability of the 
chaos-induced delocalized state as the initial state il)i{t = 
0) of the optimization procedure. Note that the target 
state is a localized Gaussian wavepacket in the left-well 
as well as the last section. Hereafter, we mainly use the 
polychromatically perturbed case (Af = 10, e = 1.0) as a 
typical controlled system i?o(i)j which is strongly chaotic 
as shown in Fig. 2. 




FIG. 7: Long-time optimal control for the delocalized initial 
state in the polychromatically perturbed case (M = 10, e = 
1.0). (a)Time-dependence of the squared final overlap Jo{t) 
after 50 iterations, (b) Convergence of the optimized objective 
functional values J and Jo as a function of the number of 
iteration. The panels (c) and (d) show the optimized fields 
E{t) and the power spectrum I{u)) of the field. We used 
Of = 1.0 and tf = 300 in the calculation. As an initial state we 
used the initial delocalized state obtained by time-evolution 
up to f = 5.54 X 10^ in the system with M = 10, e = 1.0. 

Figure 7 shows the relatively long-time control of the 
delocalized state in the polychromatically perturbed case 
(M = 10, e = 1.0). The achievement of the optimization 
is almost same as the case of the Gaussian wavepacket. 
Figure 8 shows the short-time control for the same case. 
It is difficult to design the optimal field for the short 
target-time [tf — 30). As seen in Fig. 8(b), it seems that 
the objective functional value saturates around Jo 0.8 
even for 200 iterations, and the state is trapped in a local 
minimum in the optimization process. The optimized 
field consists of a lot of complicated modes as seen in 
Fig.8(d). 

Figure 9(a) and (c) show the final states by using the 
optimized fields in Fig. 8(c) and Fig. 7(c), respectively. 
Figure 9(b) and (d) show the time-dependence of ApAg 
during the time-evolution by the optimized field. Al- 
though the quantum fluctuation gradually decreases from 
the initial large value of the delocalized state toward 
small value (~ 0.5) in the final Gaussian wavepacket, it 
is hard to squeeze the delocalized state in the short-time 
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FIG. 8: Short-time optimal control for the delocalized initial 
state in the polychromatically perturbed case (M = 10, e = 
1.0). (a)Time-dependence of the squared final overlap Jo{t) 
after 200 iterations, (b) Convergence of the optimized objec- 
tive functional values J and Jo as a function of the num- 
ber of iteration. The panels (c) and (d) show the optimized 
field E{t) and the power spectrum I{uj) of the field. We used 
a = 0.1 and tf = 30 in the calculation. As an initial state we 
used the initial delocalized state obtained time-evolution up 
to t = 5.54 X 10^ in the system with M = 10, e = 1.0. 



. (a) 


final distibution 

r 


1 
/ 




6 -4 


-2 2 4 


1 


'i final distibution 
i target state 


1 


I. 




50 100 150 200 250 300 



FIG. 9: (a), (c) The probability density of the target 
wavepacket and final wavepacket achieved by using the op- 
timal field in Fig. 8(c) and Fig. 7(c), respectively, (b), (d) The 
time-dependence of the uncertainty product ApAq during the 
time-evolution by the optimal field in Fig. 8(c) and Fig. 7(c), 
respectively. 



control. 

Figure 10 shows the final overlap Jo as a function of 
target time tj. The time-dependence of the functional 
values Jq show irregular increase with time. It seems 
that the Jq suddenly increases around tf = 100 and al- 
most saturates to Jo 1.0 around ~ 1000. Apparently 
the longer the target time becomes the better optimiza- 
tion achieves. We can expect that the optimized field 
converges to the global minimum for the longer target 
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t, 

FIG. 10: t / —dependence of the optimal control of the delocal- 
ized state in the unperturbed case. The squared final overlap 
Jo{tf) after 50 iterations as a function of target time tf. We 
used Q — 1.0 in the all calculations. As an initial state we 
used the initial delocalized state obtained by time-evolution 
up to f = 5.54 X 10^ in the system with M = 10, e = 1.0. 



time. 

It should be emphasized that even for the long-time 
control the difference of the controllability between the 
delocalized state and the Gaussian wavepacket is the 
number of iterations for the convergence that represents a 
cost for the optimization. The chaos-induced delocalized 
state consists of a lot of eigenstates on the unperturbed 
basis due to the mixing of multilevel caused by the per- 
turbation. 



VI. LOCAL OPTIMAL FIELD 

So far, we used a typical sinusoidal field = 
sini as an initial guess at the first iteration step. The 
final optimized field E{t) does not strongly depend on 
the trial field for long-time control. In this section, we 
investigate the influence of the initial trial field E^"°'^ (t) 
on the optimal field E{t) for the short-time control as 
in Fig. 8. For that purpose, as the initial state we use 
the delocalized state in the polychromatically perturbed 
double-well potential [M — 10, e = 1.0). All numerical 
settings are same to the last section except for the first 
guess of the iteration. 

In Fig. 11 we show several optimal fields for the short- 
time control in the transition between the delocalized 
state and the Gaussian wavepacket, where we used the 
following trial fields as the first guess, £2"°"' (t) = sin 2.8t, 
E'r^'^t) = ^J2fi^smin,t), El"-^it) = fit). The 
f(t) is randomly and uniformly distributed in the range 
[—1,1]. As shown in the right panels of Fig. 11 the shape 
of the power spectra corresponding to the these optimal 
fields are quite different with respect to one another. The 
function form of the optimal field E{t) strongly depends 
on the trial function form (i) as the first guess. Fig- 

ure 12 shows the convergence to the optimal fields in the 
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FIG. 11: The trial field E'"'''{t) dependence of tiie short- 
time control between the delocalized state and the Gaussian 
wavepacket in the polychromatically perturbed case {M = 
10, e — 1.0). The left panels (a,c,e) are the optimized fields 
after 200 iterations for three different trial fields as the first 
guess, Ei"''\t), Ei"°'\t) and given in the text. The 

right panels (b,d,f) are the power spectra corresponding 
to the optimal fields in the left panels. We used a = 0.1 and 
tf — 30 in the calculation. As an initial state we used the 
initial delocalized state obtained (ug by) time-evolution up to 
t = 5.54 X 10^ in the system with M = 10, e = 1.0. The final 
overlaps are less than 0.8 in the all cases. 
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FIG. 12: Comparison with the results for the differ- 
ent first guesses E^""'\t) in Fig.il. (a)Time-dependence 
of the squared final overlap Jo(t) after 200 iterations, 
(b) Convergence of the optimized objective functional values 
J and Jo as a function of the number of iteration. The other 
parameters are the same as a case in Fig. 8. The curves de- 
noted by the figures (2, 3, 4) express the results corresponding 
to the E^"°'\t), El^^^t) and ET'''\t). 



cases of the different initial fields E^"''\t). The opti- 
mized values Jo after 200 iterations are about 0.8 in all 
cases, and the control is incomplete. 

As a result, it is conjectured that during the optimiza- 



tion process the numerous local minima prevent the con- 
vergence to the global minimum corresponding to the 
target state. This difficulty in finding the global mini- 
mum comes from a feature of the chaos-induced delocal- 
ized state with the complicated phase structure like an 
entangled state, ft is expected that the transition be- 
tween the delocalized state and the Gaussian wavepacket 
can be achieved by the more complicated control due to 
existence of a lot of local minima in the solution space 
via the excited states with high energy than that between 
the Gaussian wavepackets. In the short-time control, for 
the sake of the manipulation of the control field for the 
delocalized states, we need some strategies in order to 
throw off the trap in the local minimum, such as genetic 
algorithm [s^l and chaotic itinerancy 'SS] . Further inves- 
tigation is necessary to find the global optimal field for 
short-time control problem concerning the chaos-induced 
delocalized state |4g. The development of new optimal 
control method is out of scope in this paper. 



VII. CHAOS-ASSISTED CONTROL 

In this section, we consider infiuence of chaos in the 
controlled system Hq (t) on the controllability for the de- 
localized state. All numerical settings are same to the 
last section except for the first guess of the iteration ex- 
cept for value of e in A{t). 

As shown in Sect. II the chaotic behavior in the per- 
turbed double- well system is enhanced by the increase of 
the perturbation strength e. 
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FIG. 13: e— dependence of the short-time control in the 
polychromatically perturbed case (M = 10). (a)Time- 
dependence of the squared final overlap Jo{t) after 200 iter- 
ations, (b) Convergence of the optimized objective functional 
value J and Jo as a function of the number of iteration. We 
used a — 0.1 and = 30 in the calculation. As an initial 
state we used the delocalized state obtained time-evolution 
up to f = 5.54 X 10^ in the system with M = 10, e = 1.0. 



Figure 13 shows the short-time control for various val- 
ues of e(= 0.1,0.4,1.0) in the polychromatically per- 
turbed case {M = 10). It found that the controllability 
is suppressed and the final overlap less converges to the 
target state as the decrease of the perturbation strength. 
We can regard that the strong chaotic mixing in the 
double-well system assists the control of the transition 
between the quantum states. We call the phenomenon 
chaos-assisted control. 
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A simple scenario for chaos-assisted control is at- 
tributed to the following fact. In a weak chaotic system 
with small perturbation strength, the existence of the 
dynamical structure in the phase space can work as the 
obstacles for the motion of the wavepacket. On the other 
hand, the strong perturbation can remove the restriction 
of the motion and make displacement of the wavepacket 
smoother than the weakly perturbed case. The chaos- 
assisted control can be also observed in short-time con- 
trol between Gaussian wavepackets discussed in the last 
section. 

What cases are uncontrollable even for long-target- 
time tfl The energy difference between the initial and 
target states may make difficult to control the transition 
of the delocalized state if the initial state is at least a non- 
stationary state with high energy. We have confirmed the 
difficulty of the control even for a very long-target time 
optimization when we use a delocalized state with the 
higher energy as an initial state. Then, it can be ex- 
pected that the harder the control becomes the more ef- 
ficient the chaos-assisted control may work, although the 
control is still far from the perfect one. The details will 
be reported elsewhere in focusing on the more systematic 
investigation for this point (461] . 

VIII. SUMMARY AND DISCUSSION 

In summary, we numerically investigated controllabil- 
ity of the wave packet dynamics in the polychromatically 
perturbed one-dimensional double-well system based on 
the optimal control theory. The controllability depends 
on the target time tf and the initial state. The results 
are summarized as follows. 

(1) In the cases of the relatively long-time control 
for the displacement between the localized Gaussian 
wavepackets in the wells, it can be well-controllable re- 
gardless of the perturbation type. 

(2) In the cases of the relatively short-time control for 
the Gaussian wavepackets, we can obtain the almost op- 
timal field although the convergence of the iteration is 
slowly. 

(3) When we use a delocalized state as the initial state, 
which is obtained by the time-evolution under the poly- 
chromatic perturbation {M = 10, e = 1.0), the quantum 
control becomes more insufficient than the case of the 
Gaussian wavepacket. The convergence of the iteration 
is slow because of numerous local minima. 

(4) As the increase of the perturbation strength e in 
the controlled system the controllability of the delocal- 
ized state can be enhanced because the restriction caused 
by the dynamical structure could be broken in the fully 



chaotic state.. We called the phenomenon chaos-assisted 
control. 

It is interesting that whether or not the short-time 
optimization of the delocalized state is well-controllable 
by the other control methods, for example, noniterative 
method and so on, proposed by some groups 0, 0| . 

We modeled the other degrees of freedom coupled with 
the double- well system as external coherent perturbation 
in Eq.(l). A control problem for a quantum system in- 
teracting with a bath composed of finite (or infinite) de- 
grees of freedom is interesting ,22^ ■ The application 
of OGT to laser driven hydrogen tunneling by using ul- 
trafast laser pulses have been proposed by some groups 
[36]. Then decoherence through coupling to intermolecu- 
lar and intramolecular degrees of freedom is inevitable in 
the realistic system. The influence of the other degrees 
of freedom is modeled as coupling with an infinite collec- 
tion of harmonic oscillators characterized by a broad and 
continuous distribution of frequencies. From the point of 
view, we modeled the effect of coupling with the other 
degrees of freedom as a coupling with finite collection 
of highly excited harmonic oscillators with some discrete 
special frequencies. The control of the dissipative quan- 
tum state and entangled state are directly related to the 
realization of the quantum computation 3X, ,38] . 

The transfer mechanism of the well-localized Gaus- 
sian wavepacket between the wells has an analogy with 
Mott's variable range hopping (VRH) between the lo- 
calized states in (ig a) disordered system containing the 
roughness in the potential energy landscape, which is a 
dominant transport mechanism in low-temperature limit 
[39j ]. On the other hand, the conductive prop erties of 
DNA have recently attracted interest [ifillillliali^ . The 
lambda phase DNA (A— DNA) has been sometimes mod- 
eled as a one-dimensional disordered system, and there 
has been an explanation for the electrical conductivity 
along the DNA double helix, by using the VRH based 
on the temperature dependence of the conductivity |44| . 
Therefore, the control of the electron transfer in the dis- 
ordered system composed of fluctuating multi-well po- 
tential is also interesting problem. 

Moreover, we can realize the another delocalized state 
in the polychromatically perturbed one-dimensional dis- 
ordered system "4^. The more details of the controlla- 
bility of the delocalized states in the disordered system 
and the double- well system will be reported in the future 
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